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Hadrons are composite objects made of quarks and gluons, and during a collision one can have 
several elementary interactions between the constituents. These elementary interactions, using an 
appropriate theoretical framework, can be related to the total and elastic cross sections. At high cm. 
energy it also becomes possible to identify experimentally a high p± subset of the parton interactions 
and to study their multiplicity distribution. Predictions of the multiple interactions rates are difficult 
because in principle one needs to have a knowledge of the correlated Parton Distribution Functions 
that describe the probability to find simultaneously different partons in different elements of phase 

", space. In this work we address this question and suggest a method to describe effectively the 

f*^ ■ fluctuations in the instantaneous configuration of a colliding hadron. This problem is intimately 

' related to the origin of the inelastic diffractive processes. We present a new method to include the 

. diffractive cross section in an eikonal formalism that is equivalent to a multi-channel eikonal. We 

■ compare with data and present an extrapolation to higher energy. 
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^ r-| I. INTRODUCTION 

O ■ The evolution with center of mass energy of the total and elastic cross sections in hadron-hadron collisions, and the 
(D ' properties of multi-particle production in these interactions remain an important open problem in particle physics. 
This problem is clearly of significant intrinsic interest, but it has also important phenomenological implications: on one 
hand the estimate of the properties of hadronic interactions is obviously important at LHC to model the background 
, in the search for Higgs particles and possible forms of "new physics", on the other hand the detailed properties of 
^ • hadronic interactions properties is important in cosmic ray studies. The highest energy cosmic ray particles can 
lO \ only be detected indirectly observing the extensive air showers that they produce in the Earth's atmosphere. The 
■ development of air showers is determined by the hadronic cross sections and the properties of particle production in 
hadronic interactions, and therefore the interpretation of the available (and future) data depends also on theoretical 
assumptions about hadronic interactions. 

Hadrons are composite objects containing quarks and gluons, and it is natural to relate the cross sections for 
CJP hadron collisions to the more elementary interactions between their parton constituents. The precise method to do 
this remains however an unsolved problem. 

A possible approach is provided by the so called "minijet" models where the total and elastic cross sections are 
obtained using an eikonal formalism, in terms of the quantity (n(b, s)) that has the physical meaning of the average 
number of elementary interactions at impact parameter b and cm. energy t/s (for earlier work on the importance of 
minjets see @, Q)- An attractive feature of this approach is that it allows to compute the distribution of the number 
of elementary interactions that happen in a single hadron collision. This distribution can then be used, with the 
inclusion of a few additional assumptions, in Montecarlo codes to predict properties of particle production that can 
be tested experimentally. For example it is simple to see that events with a large number of parton scatterings must 
have a larger multiplicity and a more complex structure. 

The original version of the minijet approach and several subsequent ones did not include in a consistent way in 
their formalism the inelastic diffractive processes. Following an approach introduced long ago by Good and Walker 
[J] , several authors 0, @, 0, HI, HI have indicated possible methods to include diffraction using a multi-channel eikonal 
formalism. In this work we rediscuss this problem, and suggest an alternative method to include the diffractive cross 
section in the eikonal formalism. This method is mathematically equivalent to the multi-channel eikonal method but 
offers additional physical insight. 

The fundamental physical idea to explain the existence of inelastic diffraction introduced by Good and Walker Q| is 
to assume that inelastic diffraction emerges because an interacting hadron can be seen as a superposition of different 
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states that undergo unequal absorptions. It is natural, as originally proposed by Miettinen and Punrplin [ic| . to 
identify these "transmission eigenstates" as different "configurations" of the parton constituents inside a hadron. In 
this theoretical framework the estimate of inelastic diffraction requires some understanding of the ensemble of such 
parton configurations. This appears as a daunting task. 

A possible approach is to make the dramatic approximation of reducing the space of parton configurations to a 
finite dimensional space (in fact spanned by as few as two or three base vectors) and to construct explicitly a matrix 
transition operator. An explicit example (originally constructed in [H[) of this approach is also discussed in this work. 

The alternative method we propose here is to construct a mapping from the space of the hadron configurations to 
the real positive numbers, so that the number of elementary interactions for the configuration C is n(C) = (n) a(C) 
(with (n) the average over all configurations). The probability distribution of a (together with a model for (n)) is 
then sufficient to compute the total cross section and its different components (elastic, absorption and diffraction). 

The consistent inclusion of diffraction in the theory is very important, because it changes dramatically the relation 
between the inclusive parton cross sections and the directly observable hadron cross sections. 

This work is organized as follows. In the next section we start discussing the problem of multiple interactions in a 
narrower but better defined sense, that is limiting our considerations to hard high p± interactions that on one hand 
can be easily identified experimentally, and on the other hand have an inclusive cross section that is calculable in 
perturbative QCD using the standard Parton Distribution Functions (PDF's). We will show that the calculation of 
the multiplicity distribution of these hard interactions requires additional theoretical constructs. 

In section IIIII after rewiewing some general formalism about total and elastic cross sections we discuss a standard, 
single channel, version of the eikonal formalism, that has been used in the original minijet model of Durand and Pi 
and in many other works. 

In section IIV1 after a brief introduction to inelastic diffraction we recall the basic ideas of the Good and Walker 
[J] method, and then we discuss our implementation of this multi-component ansatz in terms of the function p(a) 
(where a is a real positive variable). We call the function p(a) the "effective configuration probability distribution". 

In section [V] we present a simple form for the function p(a) that depends on a single parameter, and use this form, 
together with a (2 parameters) parametrization of the function (n(b, s)} to describe the available data on pp and pp 
scattering at collider energies. 

In section IVll we discuss the energy dependence of the parameters of our model and discuss extrapolations to higher 
energy (LHC and ultra high energy cosmic rays). Section [VIII offers some final considerations. 



II. MULTIPLE INTERACTIONS 



The problem of multiple interactions in a hadron-hadron collisions is usually discussed in the context of a calculation 
of the total (or inelastic) cross section and referring to the total number of elementary interactions in a collision. Such 
a general discussion has serious difficulties. Theoretically the concept of the "total" number of elementary interactions 
in a collision is not really well defined. For example, one usually divides the elementary interactions into two classes: 
"soft" and "hard" choosing a rather arbitrary cutoff in pj_ , however "soft" interactions cannot be considered as truly 
elementary, since they are effective processes, as for instance a "pomeron exchange" that, after decades of efforts, 
remains a somewhat elusive concept; and the theoretical "counting" of the number of interactions in a collision has 
significant ambiguities. On the other hand, experimentally it is essentially impossible to measure the number of soft 
interactions in one collision, and one can at best obtain only some partial and indirect information from the study of 
particle multiplicities. 

To avoid these difficulties, in this section we will discuss a more limited but much better defined problem, namely the 
production of parton-parton scatterings with a transverse momentum larger than a chosen threshold. If the threshold 
p™ ln , is sufficiently large (above a few Gev) these scatterings are, at least in principle, experimentally identifiable as 
pairs of back to back jets. Moreover, given the colliding hadron PDF's, the inclusive differential cross section for the 
production of pairs of jets can be estimated in perturbation theory from the well known expression: 



dp± dxi dx2 
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In this expression /j 11 (x, p?) [f£ 2 (x, /x 2 )] is the Parton Distribution Function (PDF) for parton of type j (fc) in the 
hadron hi (/12) ar the scale /j 2 ~ p\, and dajk^j'k' /dp±(p±,s) is the differential cross section for the parton-parton 
scattering of type j + k — > j' + k' at the squared cm. energy for the parton-parton interaction s = sxi %i. The 
fractional momenta x\^ are connected to the rapidities of the observed jets by the relation: 
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Integrating equation (JTJ) over the phase space region p± > p" lm and all allowed jet rapidities one obtains the inclusive 
jet cross section: 

* jct (pT\V~s) = I ^ dp ± f d Xl f dx 2 \ V f^( Xl ,fx 2 ) ft(x2,(S) deTj ^ j ' k ' (p±,s) \ (3) 

JpT" Jip'l/s Jip 2 ± Ks Xl ) (j.r~'jr dp± J 

More in general, one could consider the production of jet pairs in a more limited region of phase space selecting for 
example only jets in a certain range of rapidities, with appropriate choice of the limits in the integration over phase 
space. The discussion below remains valid also in this case, and in the following we will denote with <7j e t the inclusive 
production of jets in a fixed kinematical region determined by p™ m and appropriate cuts in the jet rapidities, leaving 
the dependence on these kinematical cuts implicit. 

The quantity erj t(s) is a cross section for parton-parton interactions, and therefore must be interpreted in the 
appropriate way. Its physical meaning is to give the inclusive cross section for the production of jet-pairs in the 
chosen kinematical region. This means that when a detector collects the integrated luminosity £i n t, the expected 
number of jet pairs is Lj nt Oj c t. Since the total number of inelastic scattering events is L lrA a- mc \, the ratio <7j e t/cinei 
is the average number of jet pairs produced in one inelastic interaction. In principle it is possible, and in fact it will 
happen for p" lln sufficiently small and/or sufficiently large, that Oj C t exceeds <7i ne i. This simply implies that some 
events must contain more than one parton-parton interaction. 

In general, an inelastic event can have 0,1,2 or more hard interactions. A natural problem is the estimate of the 
relative frequencies of events that have jet multiplicity fc. The probability pu that an inelastic event contains exactly 
fc pairs of jets can be expressed as the ratio: 

Pk = 2*- (4) 



The partial cross sections <7^ et must satisfy the sum rules 



^ (jf * = (Tinel (5) 



k = CT jet (6) 

fc 

and therefore 

(*>=£*pjr=-?* ( 7) 
k crinci 

It is important to stress that the set of partial cross sections cr^f*, or equivalently the probabilities pu are observable 
quantities. For large p™ m the identification of the hard interactions is experimentaly straightforward, however the 
inclusive jet cross section ujet is much smaller than a mc \ and the jet multiplicity distribution becomes "trivial" and 
only po and p\ arc different from zero. For p™ m sufficiently small the jet multiplicity distribution become broader, and 
the probability to find more than one hard scattering in a single event becomes significant; however at the same time 
the experimental identification of the hard scatterings becomes more difficult. At LHC it should however be possible 
to identify a value p™ m sufficiently small to result in an interesting multiplicity distribution of hard interactions, and 
sufficiently large to allow the measurement of such a distribution. 

Such an experimental study should be compared with a theoretical prediction. The calculation of the multiplicity 
distribution of the hard interactions is in fact a very difficult, unsolved problem, that, as we will discuss in the 
following, requires the introduction of new ideas, beyond the use of the standard PDF's, that are only sufficient for 
the calculation of the inclusive jet cross section. One may note that the obvious sum rule ([5]) actually implies that 
the calculation of the partial cross sections requires unavoidably a complete theory for the inelastic hadron-hadron 
cross section. 

As a first step toward the calculation of the set of partial cross sections crl et one can note that it is natural to expect 
that collisions at different impact parameter will result in a different number of hard interactions, with small (large) 
b corresponding to a larger (smaller) number of interactions. The average number of interactions at a fixed impact 
parameter b can be calculated as: 



(n iet (b,s,pT a )) = J d 2 h J d% P int (b- b!+b 2 ) x 

J d P± J dx x J dx 2 F^( Xl M,» 2 ) ^ i2 (x 2 ,6 2 ,/i 2 ) d&j ^f k ' ' (Px,3) ^ 
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(where we have left implicit the integration limits over p±, x\ and x<i). The expression (|SJ) differs from equation Q 
for three reasons: (i) it replaces the PDF's fj(x,p 2 ) with the impact parameter dependent PDF's F^(x,b, /i 2 ); (ii) 

it includes two additional integrations over b\ and b 2 that describe the positions in transverse space of the partons 
inside the two hadrons; (hi) the function Pi n t(b — bi + 62) is also included. 

The impact parameter dependent PDF's F^(x, b, p 2 ) describe the probability to find a parton of type j in hadron h 
with fractional longitudinal momentum x at transverse position b (with respect to the hadron cm.). These functions 
are related to the standard PDF's by the relation: 

fHx,p?) = [ d 2 bF?(x,b,» 2 ) (9) 



The function P ln t{b) describe the probability density that two partons (each in a different hadron) separated by the 
distance b in transverse space interact. The function has the normalization: 

d 2 b P int (b) = 1 (10) 

The simplest choice for Pi n t(b), is a delta function (P; n t(fr) = S 2 (b)), or more in general (but at the cost of introducing 
a new parameter) a gaussian of width a. 

Using equations ^ and (p~0|) it is simple to verify that integrating (nj et (6, s)) over all impact parameters one obtains 
the correct result: 

r d 2 b (n jet (M)) =0jet(s) (11) 

The physical interpretation of equation §E§ is very straightforward. To obtain the expected number of hard scattering 
in a hadron collision at impact parameter b one must integrate over the distribution of the partons in trasverse space 
around the center of mass of the colliding hadrons. The standard PDF's integrate over this transverse space variables, 
and are therefore insufficient for the calculation, and are only capable of giving the inclusive result. 

To proceed in the calculation one must obtain informations about the impact parameter dependent PDF's. The 
simplest hypothesis is to assume that the dependence on impact parameter of the F^ functions can be factorized: 

F 1 h (x,b,p 2 ) = f! j l (x,p 2 ) p h (b) (12) 
where the function ph(b) satisfies the normalization condition: 



d z b p h (b) = 1. (13) 

With this factorization assumption the quantity (n(b, s)) can be written as: 

(n jot (M)> = *jet(«) A{b) (14) 

where the "hadron overlap function" A(b) is: 



A(b) = J d A h J d 2 b 2 Phl (61) p h2 (6a) P int (b -b 1 + b 2 ) (15) 
and (using (jTUJ) and (fT3"|) ) satisfies the normalization condition: 

d 2 b A(b) = 1 (16) 



A reasonable first approximation for the overlap function A(b) proposed by Durand and Pi [l[ is to estimate it 
from the electromagnetic form factor of the colliding hadrons (see appendix [Aj). Another (more phenomenological) 
approach has been to assume a gaussian form, of a superposition of Gaussians. 

It should however be stressed that the factorization hypothesis (|12p has no serious motivation beyond its simplicity, 
and it is likely to be incorrect. In fact, the impact parameter PDF's can be in principle calculated from the generalized 
parton distribution functions (GPDF's) [H, E3, [3l . Work on this subject is in progress. 

The need to consider the dependence on the transverse degrees of freedom introduces a serious complication and 
uncertainty in the calculation of the partial cross sections cr J fc ct , but unfortunately, even having a good theoretical 
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control of the overlap function A(b) (or better a detailed knowledge of the impact parameter dependent PDF's) is not 
sufficient to estimate the partial cross sections. 

To complete the calculation one must make some hypothesis about the fluctuations in the number of hard inter- 
actions for collisions at a fixed impact parameter (and cm. energy). The simplest hypothesis is to assume that 
the fluctuations are simply Poissonian. The partial jet cross sections are then calculable integrating over all impact 
parameters the Poisson probability: 

4%) = J d 2 b ( " jct( fc 6 ! ' s))fc e -<«^(M> (17) 

The hypothesis that the multiplicity distribution of the hard interactions at a given impact parameter is Poissonian 
is however not necessarily correct, and some simple considerations suggest that in fact this distribution is considerably 
broader than a Poissonian. An argument in this direction can be developed as follows. At the instant t a hadron can 
be described as an ensemble of partons each having a certain longitudinal momentum x and a trasverse position b. 
The set of values {qj, aJj, bj}(j=i,jv) {Qj is the complete set of quantum numbers of each parton) is the "instantaneous 
configuration" of the hadron and will be denoted with the symbol C. The probability to find hadron h in a certain 
instantaneous configuration C can be denoted as Pfe(C) with the normalization condition: 



JdCP h (C) = l (18) 



where the integration over dC indicates formally the sum over all possible configurations. 

The interaction between (for example) two protons, in the cm. frame, lasts a crossing time i C ross, that can be 
estimated as: 

2Rp _ ARprrip 

Across — r— V ^/ 

7cm V s 

where R p is the linear size of the proton of order 0.5 fm. Because of the Lorentz length contraction the crossing time 
shrinks oc s~ x / 2 

The time required for the radiation and absorption of partons, that is for changing the parton configurations is of 
order R p and is therefore much longer than the interaction time, ft is therefore a reasonable approximation to assume 
that during one interaction the hadrons appear as "frozen" each into its own configuration. 

The expected number of hard interactions in a collision of impact parameter b is determined by the parton con- 
figurations of the two colliding hadrons. If the configuration is composed by few hard partons, as for example is the 
case when most of the hadron energy is carried by the valence quarks, the number of hard scattering is suppressed, 
while if the configurations contain many gluons with fractional energy of order x ~ (p™ m ) 2 /s the number of hard 
interactions is enhanced. 

We will denote as n- }e t{b, Ci, C2) the expected number of hard interactions in a collision at impact parameter b 
between two hadrons with configurations Ci and C2. Integrating over all possible hadron configurations one must 
recover the expected value of the jet multiplicity, therefore we can write formally: 

fdd J dC 2 Ph.id) P h2 (C 2 ) n jct (6,d,C 2 ) = (n jet (M)) • (20) 

If we make the further assumption that the actual number of interactions in a single collision at impact parameter 
b with the hadrons in the configurations Ci and C 2 has a Poisson distribution around the expected value, we can now 
write the partial jet cross sections as: 

of = J d 2 b J dCt J dC 2 P hl (d) P h2 (C 2 ) I [^jct(6,Ci,C 2 )] exp[ _ n . ct(5;Cl;C2)] | (21) 

Of course equation (|2ip is only a formal solution of our problem since we have not yet developed the instruments to 
estimate the probability of the different configurations and to perform the integration over the configurations space. 

In order to make progress we will make the assumption that the expected value of the jet multiplicity for a certain 
configuration of the hadrons is related to the averaged one by the factorized relation: 



rajet(b,<Ci,C 2 ) = (rijet(M)) «(Ci,C 2 ) 



(22) 
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where a(Ci, C2) is a real positive number. With this simplifying assumption we can construct a function p(a) that is 
independent from b: 



P (a) =JdC 1 J dC 2 P hl (d) P h2 (C 2 ) <S[a(d,C 2 ) - a] (23) 
It is straightforward to see that because of the normalization condition (|18p one has: 

da p(a) = 1 (24) 



consistently with the interpretation of p(a) as a probability density; moreover because of equation (|20p one has also 
that: 



n 



da a p{a) = 1 (25) 
Using the definition of (|23p one can rewrite (|21[) in the simpler form: 



jot - ^6 r da p(a) [ { " j * (fe ' S))k cxp [-a (n jot (b, a))] 



In this equation we have finally arrived to write the jet multiplicity distribution in terms of (nj e t(6, s)) and the function 
p(a), that is unknown, but has fixed normalization and first moment. The expression for the partial jet cross sections 
given in equation (|17[) ) that was calculated ignoring the effects of fluctuations in the parton configurations is recovered 
in the limit where the function p(a) has vanishing width and reduces to the form p(a) = 6 [a — 1]. 

It is straightforward to see (using equation (|25p) that one has: (fc) = crj et / c"inei m agreement with equation ([7]). The 
second moment of the jet multiplicity distribution can be expressed as: 

(fc 2 ) = J- £ fc 2 af = + £±™l [ dH (n iet (b, S )) 2 (27) 

^"incl 7 ^inel <^incl J 

k 

where 

w = da (a 2 - 1) p{a) (28) 



Jo 

is the variance of the p[a) distribution. Using the factorization hypothesis (|14p one can rewrite equation (|27[) as: 

( A 2) = £JSL + (1 + 1£? ) Hi*. fd 2 b[A(b)} 2 (29) 

Cine! Cinol J 

The width of the jet multiplicity distribution is therefore determined by the geometry of the hadronic matter (that 
is from the shape of the overlap function A(b)) and from the variance w of the function p(a). Note that the simple 
cikonal model corresponds to a vanishing w and therefore to the minimum possible (fc 2 ) (for a fixed overlap function 
A{b)). 

The partial cross sections oi e and the quantity (fc 2 ) are (at least in principle) measurable with observations of 
the jet multiplicity distribution. From these measurements one can obtain information about the properties of the 
function p(a). 

In fact we have already some (indirect) information about the jet multiplicity distribution, and there are indications 
that predictions based on the simple eikonal approach are not adequate. Some of the most sophisticated Montecarlo 
instruments for the study of high energy hadron collisions, such as PYTHIA [lj| [l6| and HERWIG [l?], [3 include a 
treatment of multi-parton interactions following the simple eikonal model. The algorithms for multi-parton interac- 
tions in these Montecarlo codes compute an inclusive jet cross section a- ]e t, assume an overlap function A(b), and them 
generate a number of elementary interactions according to the probability distribution of equation (|17p or equivalently 
of equation (|2l))) with p(a) = S[a — 1]. In most cases the presence of multiple interactions cannot be experimentally 
resolved but is detectable statistically, for example from fluctuations in the charged particles multiplicity distribution. 
In order to reproduce the broad fluctuations of the data, it appears necessary to construct some ad hoc functional 
form for the overlap function. For example PYTHIA [3, [l6[ uses a double-Gaussian with a denser core and more 
extended halo. The overlap function A(b) however (as will be discussed in more detail in the following) in principle 
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also determines the t distribution of elastic scattering, and therefore one does not have the freedom to modify it in 
an arbitrary way. 

The introduction of a non vanishing variance for the function p(a) allows to modify the width of the jet multiplicity 
distribution without modifications of the overlap function A(b). 

One open problem in high energy hadron collisions, that is likely to be relevant in the interpretation of future 
data on possible manifestations of new physics in high energy collisions is the problem of the so called "underlying 
event" . The observations [li^ . |20| show that in events with the presence of a high p± scattering, the "environment" 
that accompany the observed jets has an average transverse momentum higher than what is found in minimum bias 
events. Montecarlo codes like HERWIG, ISAJET or PYTHIA at present describe only partially these effects [HHll- 
The origin of a higher "ambient level" of p± for the underlying event can be related to the presence of additional (softer 
and unresolved) parton scatterings that accompany the observed jets. The correct description of the underlying event 
is therefore related to a good theoretical control of the multiplicity distribution of parton interactions. A Montecarlo 
inplemcntation of multi-parton interactions based on an impact parameter picture for hadronic collisions naturally 
contains some of the qualitative features observed in underlying events, because events selected with a high p± jets 
are more likely to be central collisions, and therefore are more likely to contain additional parton scatterings. A non 
vanishing variance of p(a) should however enhance the differences between underlying and minimum bias events, and 
could therefore play a non negligible role in the description of the data. A quantitative study of this problem with 
Montecarlo methods is a goal for future work. 

The problem of the fluctuations in the configurations of partons in a hadron is also intimately related to diffractive 
scattering, as we will illustrate in the following section. 



III. TOTAL AND ELASTIC CROSS SECTIONS 



A. General Formalism 



In the following we will need to consider also the elastic and total cross section in hadronic scattering. The elastic 
scattering amplitude in the collision of two hadrons of type hi and /12 at c.m. energy y/s can be written [2l[ as a 
(2-dimcnsional) integral over impact parameter: 

— e^r o i(M (30) 

27T 

where r e i(6, s) is the profile function, that without loss of generality can be written as: 

r cl (M) = 1 - e~ x(M (31) 

with \(b, s) the eikonal function. In the notation we are leaving implicit the dependence on the type of hadrons 
participating in the collision. The clastic scattering amplitude is related to the differential cross section by: 

^(*' s ) =7r $f(£ s ) =7r \MV=i,s)\ 3 (32) 

In this equation t = [ph x — p'ht) 2 1S the squared momentum transfer, and we have used the approximation t = — \q\ 2 
that is valid for small —t and high energy. Integrating over all t one obtains: 

<rei(«) = / J 2 b |r 0l (M)| 2 (33) 

The total cross section is related to the imaginary part of the forward elastic scattering amplitude by the optical 
theorem and is given by: 

otrt(8)=47r9[Fei(0,*)]=2 / d 2 b R [T el {b, s)} (34) 

Combining equations (|33[) and (|34[) one obtains an expression for the inelastic cross section as an integral over the 
impact parameter: 

a ine i(s) = J d 2 b {l - |1 - r e i(6, s)| 2 } (35) 
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Equation ([32]) gives the exact shape of the differential elastic cross section. For small \t\ this shape is to a good 
approximation a simple exponential [da c \/dt cx e Bt ) and it is convenient to define the slope B e \ of elastic scattering: 



Bel(s) 



/ da e \ N 




( da c \ \ 


V dt , 




\ dt J 



(36) 



ft is straightforward to sec that if the profile function is real, the slope B e \ can be calculated as: 



Bel (a) 



J d 2 b b - r el (b,s)\ x y^br 



:,<M)> 



and measures the value (b 2 ), of the profile function. 

Using the approximate exponential shape of the differential cross section one finds the relation 

|F cl (0,s)| 2 (3[F cl (0, S )]) 2 \ l + P \ 2 

a e \(s) ~ 7T ; ; = 7T 



B e l(s) 



B e l{s) 



(37) 



(38) 



(where p is the ratio of the real to imaginary parts of the forward elastic amplitude). From the optical theorem (|34[) 
one obtains the relation: 



16 7TB, 



cl 



(39) 



B. The simple eikonal model 

The expressions (|301l35| allow to compute the total and elastic cross sections for hadron collisions in terms of the 
profile function r e i(6, s) or equivalently of the eikonal function x(b, s). We remain with the task of constructing these 
functions. Physical insight on T c \(b,s) and x(6, s) can be obtained using a well known analogy with the classical 
treatment of the absorption and scattering of a plane wave of light from an opaque screen. If the ratio between the 
amplitude "just behind" the screen and the amplitude of the incident plane wave is T(b) (where b is a 2-dimensional 
vector spanning the screen), then it is simple to obtain expressions for the total, elastic and absorption cross sections 
that are formally identical to equations (|30l [35| . In particular the expression (|35[) for the inelastic cross section 
suggests to interpret the quantity 1 — |1 — T(&)| 2 as the probability to absorb the wave at the position b on the screen. 

At high energy the elastic scattering amplitude, in reasonably good approximation, is purely imaginary, and ac- 
cordingly the profile and eikonal functions are purely real. Neglecting the real part of the scattering amplitude the 
optical analogy is then sufficient to express the profile function as: 

r cl (6, a) = 1 - y/1 -P abs {b, s) (40) 

with P a bs{b, s) the absorption probability. 

In the following we will use the approximation to consider the elastic scattering amplitude as purely imaginary, and 
the profile and eikonal functions as purely real. The analiticity of the elastic scattering amplitude that is necessary 
to respect causality, can be imposed estimating the real part with a dispersion relation. 

Equation (|40|) can be re-expressed in the form: 

r cl (fe, a) = 1 - = 1 - v^(M) = 1 - exp 

In this equation Pq = 1 — P a bs = e _ ^™' 6,s ^ is the probability of no absorption in a collision at a certain b and s. The 
notation Po{b, s) = e~( n ( b ' s ^ is motivated by the physical ansatz that "absorption" in a hadron collision corresponds 
to at least one interaction between the parton constituents of the two hadrons. If (n(b, s)) is the expected number of all 
elementary interactions between partons, and if the fluctuations of this average number of interactions is Poissonian, 
one obtains equation (|4ip for the profile and eikonal function. 

The ansatz outlined above connects the profile function r e i(6, s) to {nip, s)) that is the average number of elementary 
interactions at impact parameter b and cm. energy ^/s. In the following we will call this model as the "simple eikonal 
model". This ansatz was introduced by Durand and Pi in [l|, who proposed to compute the quantity (n(b, a)) as: 



KM)) 



(41) 



(n(b 7 s))=a jct ( S ,pf n ) A(b) 



(42) 
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where <7j e t(s,p™ ln ) is the quantity discussed in the previous section, that is the inclusive cross section for the production 
of jet pairs above a certain (fixed) p±, and the energy independent geometry factor A(b) gives the overlap of hadronic 
matter. For pp interactions Durand and Pi (see discussion in appendix |A|) chose: 

b 3 „ ( b 



^ (6) =9671f A3 UJ (43) 
with R~ 2 = 0.71 GeV 2 . 

The model in this simplest form (that has in fact a single parameter, the value of p" lm ) soon proved to be inconsistent 
with the data. Perhaps the main difficulty with the original formulation of the model is that it predicts an incorrect 
relation between <xtot(s) and B e \(s). This problem in fact is of a general nature and indicates that in the simple 
cikonal model, the factorization of (n(b, s)) as the product of two functions one of s and the other of b cannot fit the 
data. 

To falsify the factorization hypothesis Q14[) it is in fact sufficient to measure both a tot and the slope of elastic 
scattering B c \ (or more in general, because of equation (|39[) and the smallncss of p, two out of the three quantities 
Ctot, Coi and B e \) at two different values of i/s. It is always possible to find a profile function A(b) and a value of the 
eikonal function <7 e ik(si) that reproduce the observations of a to t and B D \ at cm. energy -y/sl, but then one remains 
with a single parameter er ik(s2) that can be chosen to reproduce <7tot(s2) or i? i(s2), however the two values will 
coincide only if the factorization hypothesis is valid. 

Models based on an eikonal model and the factorization hypothesis (TT41 naturally predict a correlation in the 
growth of otot(s) and B c \(s) with energy, however the growth of B e \(s) in the data is faster than the predictions, as 
was realyzed very early by Durand and Pi themselves [22j. 

This problem is illustrated in fig. [1] where we show, in the plane (a t ot,B), experimental data obtained for pp 
scattering at the ISR [23[ and for pp scattering at CERN [24| and Fermilab [25|, [26|, [27|, and compare with the 
predictions of the simple eikonal model. The thick curve is the prediction of the simple eikonal model using a 
factorized form of (n(b, s)) and using for A(b) the expression of equation (|4"5)) . This line passes through the ISR data 
points, but fails at higher energy. The other lines in the figure are calculated using the same functional form for A(b) 
but replacing the parameter R p with ?* = 1.1 R p and r = 1.2 R p . It is clear that, in the simple cikonal model the 
width in impact parameter of the function (n(b, s)) cannot be energy independent and must grow with s in order to 
reproduce the data. 

An additional problem for the Durand and Pi model, beside the one we have just discussed connected to the 
relation between atot and B c i(s) 7 is also that the energy dependence of crtot(s) of the model, driven by the growth of 
(Tjet(s,p™ ln ) with s for p™ m fixed is significantly faster than the data. 

The ansatz of the simple eikonal model has been considered by several other authors [H, [29|, HH, HH , that 
have constructed different models for the function (n(b, s)), abandoning the factorization hypothesis. For example, it 
has been suggested to decompose the function (n(6, s)) in the general form: 

(n(b, s)) = (Tcik(s) A(b, s) = 0soft(s) A so { t (b, S) + <7hard(s) ^hard(&, s) (44) 

where the two terms describe a "hard" contribution that can be calculated with perturbative methods, and a "soft" 
non perturbative part. As discussed above, the width in impact parameter of the combined overlap function A(b, s) 
must increase with s to reproduce the relation between B c \(s) and <7tot(s). Several justifications of this growth have 
been offered, however in our opinion none is really convincing. 

The energy dependence of cr e ik(s), and its decomposition in a soft and a hard part have also been the object of 
considerable discussion. In this work we will not enter in this discussion, because our main purpose here is to discuss 
the limitations of the simple eikonal model. In the simple eikonal model the quantity (n(b, s)) determines completely 
the profile function, and therefore the total and elastic cross sections. We will also consider the quantity (n(b, s)) 
attributing to it the same physical meaning that it has in the simple cikonal mode, namely the average number of 
elementary interactions at impact parameter b and cm. energy y/s; however in order to describe in a consistent 
way inelastic diffraction, we will propose a different method to connect (n(6, s)) with the profile function and to the 
observable cross sections. 

In the following we will parametrize the function (n(b, s)) as: 

<»*•»-'- W{967RWM4)} <45) 

with an overlap function that has the same form as in equation (|43[) . but with an s dependence obtained substituting 
for R p an energy dependent parameter ro(s): 

The simple eikonal model ansatz outlined above does not include a treatment of the inelastic diffraction processes. 
This is a serious limitation, because the consistent description of these processes is in fact essential. This problem is 
addressed in the next section. 
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IV. INELASTIC DIFFRACTION 

Inelastic diffraction [34| produces three classes of events: beam, target and double diffraction. In beam (target) 
diffraction the beam (target) particle is excited to a higher mass state of the same quantum numbers (with the possible 
exception of spin) while the other initial state hadron remains unchanged. In double diffraction both colliding hadrons 
are excited into higher mass states. The three processes can therefore be described as: 

hi hi — > hi hi (Beam Diffraction) 

hi h 2 -> hi h* 2 (Target Diffraction) (46) 
hi hi — > h* h* 2 (Double Diffraction) 

The experimental selection of events that belongs to these three classes is not trivial, especially for the double 
diffraction ones. These events are characterized by a large "rapidity gap" that separates particles produced in the 
decay of the two excited states h\ and h\. A complete description of beam (target) diffraction is given by the double 
differential cross sections d 2 a/(dM 2 dt) (d 2 a/ (dM 2 dt)) that gives the probability to produce an event where the beam 
(target) particle is excited to a state of mass Mi > mi (M 2 > m 2 ) with transfer momentum t. The double diffraction 
cross section is described by the 3-timcs differential cross section: d 3 er/ (dMf dM 2 dt). For single diffraction one has a 
reasonably accurate picture of the M and t dependence of the cross section. The t dependence has the characteristic 
exponential behaviour of all diffractivc-like processes, for the mass dependence the cross section grows very quickly 
above the threshold (M = m p + m^ for proton excitation), oscillates following the structure of the resonances with 
the quantum numbers of the hadron in question and then falls at higher masses roughly proportionally to M~ 2 . In 
the following we will indicate as "diffraction" the sum of these three classes of inelastic diffractive events. 

We take the point of view that any parton-parton scattering implies an exchange of color, and therefore, in the 
treatment of the cross section given above, all events with one or more parton-parton interaction must be considered 
as non-diffractive. In the simple eikonal model the total cross section is therefore decomposed into an inelastic non- 
diffractive part (with at least one parton-parton scattering) and the elastic part, with no room left for diffraction. 
This lack of inclusion of inelastic diffraction is an important conceptual problem for the simpc eikonal model. 

A. Good and Walker model 

The fundamental idea for the description of inelastic diffraction in hadronic collisions has been introduced long ago 
by Good and Walker Q as an analogy with the scattering of polarized light on a bi-rcfringent absorbing medium. 
It is well known that if a plane wave of light impinges on an absorbing screen with grayness profile T(6) one has 
absorption (oc 1 — |1 — T(fe)| 2 ) and elastic scattering ( cx |T(6)| 2 ) with a very forward diffraction pattern that depends 
on the geometry of the screen. 

If the screen absorbs differently the light polarization states, in general an incident beam of a given initial polarization 
will result in absorption, and in scattered light with the same (elastic scattering) and the orthogonal (inelastic 
diffraction) polarization state. For a more explicit example, let us consider an incident beam of light in the linear 
polarization state \x), and an absorbing screen that has the grayness profile iy (6) for light with polarization \x') and 
the grayness profile T y i(b) for light with polarization \y') with 

\x ) — cos tp \x) + sin if \y) 
\y ) = — snitp \x) + cosip \y) 

It is clear that one will have an absorption cross section according 1 — |r x /(&)| 2 — |r y /(6)| 2 and after the screen one 
will find scattered waves with both the \x) and \y) polarizations. The cross section for elastic scattering and inelastic 
diffraction turn out to be: 

a x = J d 2 b 1 1 — [1 — IV (&)] cos 2 tp - [1 - IV (&)] sin 2 ip\ 2 (47) 

a y = J d 2 b cos 2 tp sin 2 tp |iy (6) - r^(6)| 2 (48) 

Note that inelastic diffraction is non vanishing only when iy ^ T y i, and if <p is not a multiple of 7r/2, that is if 
the states {\x), \y)} do not coincide with the states {|x'), \y')} that are the eigenstates of the transmission across the 
screen. 
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It is straightforward [J] to generalize this elementary example. We can consider two complete sets of orthonormal 
states {|y m )} and {IV^)}- The states \<p m ) are directly observable, with the label m describing the invariant masses 
of the two final hadronic states (in the forward and backward hemisphere) and their particle content. Without loss of 
generality we can assume that the state \tpi) with m = 1 corresponds to the initial state of the scattering (for example 
in case of a irp collisions to the state 1 7rp) ) . 

The states are eigenstates of the scattering matrix. That is, defining as usual the S matrix as S = I + i T one 
has: 



The relation between the two orthonormal bases is given by: 

\Vm) = ^2 C'mj \4>j 
3 



(49) 

(50) 
(51) 



It is now possible to write the different components of the cross sections as integrals over the impact parameter 
dependence of the tj(b). The absorption cross section is: 



Cabs = J d 2 b 

The cross section for diffraction into the state m is: 



(52) 



d 2 b 



Y,C* m] C l3 t 3 (b) 



(53) 



The elastic cross section corresponds to a\ and therefore is 



a c \= Id b 



£|Cy| 2 *#) 

3 

Summing over all states m and using the orthonormality of the states one obtains: 



(54) 



(55) 



The inelastic diffraction cross section can be obtained subtracting equation (|5^|) from (|55|) . The total cross section 
can be obtained summing the elastic, diffractive and absorption cross sections obtaining the result: 



(56) 



in agreement also with the optical theorem. 

The problem with the approach that we have just outlined, is that it is only formal, and requires the introduction 
of many, in fact infinite parameters; moreover the nature of the eigenstates \ipj) is not physically obvious. 



B. Partonic Interpretation of the scattering eigenstates 

Miettinen and Pumplin [l(| have made the proposal to interpret the states \tf>j) (the eigenstates of the T matrix) 
as parton co nfig uration states. In their paper they also include a simple explicit methods to construct these partonic 
models (see [35l 36j for a recent rediscussion) . 

In this work we will develop the idea that one needs to integrate over the parton configurations of the interacting 
hadron. The discussion that we have outlined in section |TT] can also be applied to the current problem. Using the 
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notations that we have introduced in section pi] the label j of the T matrix eigenstates \tpj) corresponds to the direct 
product of the configurations Ci and C2 of the colliding hadrons, and the squared coefficient | Ci j | 2 corresponds 
to the probability P^ 1 (Ci) x -Ph 2 (C 2 ) to find the two hadrons in a certain configuration. One then has also the 
correspondence: 



dCi / dC 2 PhdCi) A 2 (C 2 ) 



The transmission eigenvalues tj(b) have the partonic interpretation as: 



tj (6) = 1 — exp 



nj(b) 



= 1 — exp 



n(6,Ci,C 2 ) 



(57) 



where as in section |TT] n(b, Ci, C2) is the expected number of interactions among partons in a collision with impact 
parameter b when the colliding hadrons are in configuration Ci and C 2 (the dependence on the cm. energy has been 
left implicit). The difference with respect to the discussion made above is that in the previous case one was discussing 
only a subclass of (hard) interactions, while here one refers to the total number of parton interactions. 

At this point we can again make the factorization hypothesis, that the fluctuations in the number of interactions 
at different impact parameters, being related to the distribution of parton configurations is independent from b, and 
therefore one has: 



dCi / dC 2 P hl (d) P/> 2 (C 2 ) exp 



n(b, Ci,C 2 



da p{a) exp 



(n(b))a 



(58) 



As before the function p(a) satisfies the two integral relations: 



f 

Jo 



da p{a) 



da a p{a) 



'0 Jo 
Using the identities (|57|) and (]58|) we can now rewrite in manageable form equations (|52I54I55I56II as: 



d 2 <J a bs 
d 2 b 



= 1 - 



da P (a) e -MM>" 



(59) 



d 2 a c i 
d 2 b 



da p{a) ( 



1 - e 



("(■>, 3)) < 



(60) 



d 2 (Tdiff+cl 

d?b 



da p(a) 



1-e 



(61) 



d 2 b 



= 2 



da p{a) I 1 — e 



WMli 



The elastic scattering amplitude is given by: 



d 2 b 



F el (q,s) = i / — S- h \ dap{a) 



2vr 



1 - e 



<n(b, 3 )> 



(62) 



(63) 



It is important to note that in the limit p{a) — > S[a — 1], that is in the limit where one neglects the effects 
of different parton configurations, one has that equations (f6"0"|) and (|6~lj) , that describe the elastic, and elastic + 
diffractive cross sections become identical, that is inelastic diffraction vanishes. Moreover, in this case, the expressions 
for the absorption, elastic and total cross section coincide with the expressions of the simple eikonal model that 
neglects inelastic diffraction. 

Equation (|61|) sums over all diffractive channels, and therefore loses all information about the distributions of the 
excited masses. It remains however possible to compute the t distribution for elastic scattering and for inelastic 
diffractive scattering (summing over all possible open channels): 



dOdiff+d 

dt 



Eda m 
~dT 



dMi / dM 2 



d : '' 



Odiff 



dMi dM 2 dt 



(64) 
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The differential cross section for clastic scattering da c \/dt can be calculated from equations (|32|) and (|63)) 

l 2 



dt 



db b J (by/\t\) / dap(a) ( 



1 - e" 



(n(6,s)> 



Similarly, the differential cross section for diffraction plus elastic scattering can be calculated as: 



rfgdiff+cl 

dt 



7r / e?a p(a) 



db b J (b y/\T\) (l- e - M ^ ) 



(65) 



(66) 



It is straightforward to see that for p(a) = <J[a — 1] expressions (|65|) and (f66]) become identical and equal to the 
well known expression for the simple eikonal model: 



da c \ 
~dT 



simple 



dbb J (&Vl*l) \ l - e 



(67) 



The slopes at |i| = (£? c i and Bdiff) can be calculated from the definition (j3~7| (and the analogous for inelastic 
diffraction). For example, using: 



lim — 

*->o dt 



MW\t\) 



one obtains for B P 



B r . 



db 



b 3 



da p 



(«) ( 



1- e" 



4 



da p(a) ( 1 — e 



(»M) 



(68) 



(69) 



A similar expression for -Bdiff is easily derived. 



C. Comparison with multi— channel eikonal formalism 

A method to apply to Good Walker ansatz is to construct explicitly a a transition operator as a matrix of n x n— 
dimensions. In this way the formal indices m and j in section [IV Al become simply integer indices running from 1 to n. 
To show the mathematical equivalence between this multi-channel method and the use of the effective configuration 
probability distribution p{a) is straightforward. In the multichannel approach the profile function r c i(6, s) is replaced 
by the (n x n) matrix T(b, s) that can be expressed in terms of the eikonal matrix x(b, s): 

f(M) = l-exp[-x(M)] (70) 
The eigenvalues of F (r^) and x {Xj) can be written in the form: 

MM) 



= 1 - e Xj = 1 - exp 



(71) 



The interpretation of Tj as the effect of absorption suggests that (n(b,s)} and all a^'s are real and positive. In the 
optical analogy one interprets the quantity: 

1 - [1 - r/ = 1 - e" 2 ^ = 1 - e -WM)) *i 

as the absorption probability for the eigenstate \ tpj) that corresponds to eigenvalue ctj. Assuming a Poisson probability 
distribution, the quantity [(n(b, s)) aj] can then be interpreted as the average number of elementary interactions for 
the eigenstate \ipj). Without loss of generality, reabsorbing a constant in the definition of (n(b,s)}, one can impose 
the constraint 

j 

where the quantities Pj's: 

Pi = I (V^ I ^initial) | 2 (73) 
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measure the probability overlaps between the initial state |(/?initiai) and the eigenstates \tjjj) of the transition matrix. 

The normalization condition (|72p allows to interpret the quantity (n(b, s)) as the average number of elementary 
interactions (at impact parameter b and cm. energy y/s) for the initial state | (^initial)- It is then possible to define 
the function p(a) as 

3 

and to use this function to express the results for the total, elastic, absorption and diffractive cross sections as integrals 
over a according to equations (f5T)H6"2"l) . 

It can be instructive to consider an explicit model that implements the multi-channel model. The minimum model 
that includes the 4 distinct processes of elastic scattering together with target, projectile and double diffraction, 
must obviously consider a 4-dimensional Hilbert space spanned by 4 physical eigenstates \<p m ) (that without loss of 
generality can symbolically be labeled as: \trp), \irA), \pp) and |pA)). 

The structure of this most general Good- Walker model with 4-channels has been presented in [ll[ and is discussed 
in detail in appendix [B] The most general 4-channel eikonal (for a fixed y/s) is described by the impact parameter 
multiplicity distribution (n(b, s)) and by the 4x4 matrix M that in the most general case (where the two initial state 
hadrons are not identical as in irp scattering) is defined by 4 real parameters (in the case where the two interacting 
hadrons are identical, as in pp scattering, the parameters reduce to 2). 

The eigenvalues and eigenvectors of the most general matrix M are easily calculated (as show in appendix |B|) . 
Having solved this diagonalization problem, it is simple to obtain the cross sections for the total, elastic, absorption 
and diffractive scattering. The cross sections depend on (n(b, s)} and on the 4 (or 2 for identical initial particles) 
parameters of the matrix M. The cross sections can be recast in the form (|59H62[) as integrals over a defining p(a) 
according to equation (|74[) (where the summations over j now runs from 1 to 4). 

Note that in a n-channcl eikonal, the inelastic diffractive cross sections is obtained as an explicit sum of n — 1 
terms: 

Odiff = Gm (75) 

(we are identifying the state m = 1 with the initial state) . In a realistic discussion the index m should run continuously 
over all possible diffractivcly excited states. For the 4-channel model the 3 states can be identified as representing 
target, projectile and double diffraction. 

It is interesting to study, in the framework of this most general 4-channel model, the relative importance of single 
versus double diffraction. Considering for simplicity the case of pp collision, the 4-channcl Good Walker model has 
two free parameters {(3 P and e p ). The relative importance of the elastic and diffractive processes can vary significantly 
with variation of the model parameters, and similarly the ratio between the double and single diffraction cross section 
can assume different values, however numerical studies show that to a good approximation (for the scattering of 
identical particles) one has: 

0td °bd odd odd mr\ 
= ~ = (76) 

Ocl Ccl otd ct bd 

that is the ratios of cross sections single (beam or target) diffraction to the elastic one is approximately equal to the 
ratio of the double diffraction cross section to the single diffraction one. This result can also be understood from 
inspection of the structure of the matrix M . For negliglible e p for example, the relative importance of the cross 
sections for elastic, single diffractions and double diffraction scattering are in the ratio 1 : P : 0i. This result is 
compatible with the available data on double diffraction (taking into account the large errors). 

The separate calculation of the different components of the diffractive cross sections is a significant merit for the 
4-channel model. A limit of such an approach is that it predicts the multiplicity distribution of the elementary 
interactions as the superposition of 3 (for the scattering of identical particles) or 4 (in the general case) Poissonian 
distributions of different average values, such distribution might be not sufficiently smooth for a realistic comparison 
with data. The approach of using the function of a real positive variable p(a) allows to consider implicitly an infinity 
of inelastic channels. 

In the next section we will propose a simple paramctrization of p{ot) that depends on a single parameter. 

V. EXPLICIT MODEL 

Equations (T55H6T|) allow to compute the total, elastic, diffractive and absorption cross sections in terms of the impact 
parameter multiplicity distribution (n(b, s)) (the average number of elementary interactions at impact parameter b 
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and cm. energy s/s) and of the function p(a). Unfortunately the shape of the function p(a) is not determined, all we 
have been able to establish is that the first two moments (a k ) with k = 0, 1 must be unity. It is however reasonable 
to expect that the most important property of the function p(a) is its second moment (a 2 ) or equivalcntly its width 



Given this lack of knowledge about the shape of p(a) we have chosen for it a simple analytic expression, that allows 
easy manipulations: 



p(a) 



1 



exp 



(77) 



The k th moment of the distribution is: 



therefore one finds: 



da a k p(a) 



™ fc r(fc + i) 
r (£) 



(78) 



a u = 1, 



(a) = 1, 



= 1 + UK 



(79) 



Therefore the parameter w describes the variance of the a distribution. The integer values moments for n > 3 can be 
written also as: 



(a n ) = (l + ttf)...(l + (n-l)u>) 



(80) 



Numerical examples of the function p(pt) are shown in figure [2] 

An attractive property of the functional form (|77[) for p(a) is that it allows to perform analytically the integrations 
over a in equations (|59M62p to obtain the quantities of interest. The profile function r e i(6, s) becomes: 



r c i(&, s) 



da p{a) 



1 - e" 



MM) • 



= 1-1 



(n(b, s)) w 



(81) 



The expressions for the total, elastic and diffractive cross sections become: 



d 2 otot 
d 2 b 



= 2-2 1 



»(M)>' 



(82) 



d 2 <7 c \ _ ( f 1: (n(b, s)) w 



(83) 



^ = (1 + (n(M)) - (l + M^)^ (84) 

Using for (n(6, s)) the parametrization of equation (jl5| the expressions ([52115^]) allow to compute the different 
components of hadron-hadron interactions for any given value of y/s in terms of three parameters: (c e ik,ro,w). It 
should be noted that the simple eikonal model corresponds to the case w — > and is therefore included as a limiting 
case of our model. 

As a critical remark we note that the qualitative idea behind our one-parameter modeling of p{a) is that the most 
important feature of p{a) is its second moment (a 2 ) = w + 1. This however is only true in first approximation. 
Functions p(a) that differ only for moments (a k ) with k > 2 can also produce different cross sections. An example 
of this behaviour is the 4-channel model of appendix [Bj For pp interactions this model has two free parameters (f3 p 
and e p ). The variance of p(a) in the model is (1 + f3 p ) 2 — 1 is uniquely determined by (3 P , however the value of the 
different cross sectios (elastic, diffractive and absorption) depend on both of the model parameters. 

The single parameter description of p(a) of equation (|77[) seems in any case a reasonable form to investigate 
phenomenologically the consequences of equations 
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A. 



Parameter determination 



Equations (f82H84|) allow to determine (at a certain y/s) the set (cr e ikj r o, w) of the three parameters in the model 
from measurements of (<7 tot , B, cr d iff). As an example of this parameter determination we discuss here in some detail 
the measurements performed at one particular value of the cm. energy (y/s = 1.8 TeV) by one detector (CDF at the 
Fermilab J>p collider). The CDF experiment [2{||27| has measured cr to t (l + j° 2 ) = 81.83±2.29 mbarn (that estimating 
p ~ 0.15 corresponds to cr tot = 80.03 ± 2.24 mbarn); an elastic cross section a c \ = 19.7 ± 0.85 mbarn. and a slope of 
the forward elastic cross section B c \ = 16.98 ±0.25 GeV~ 2 . In addition the CDF collaboration has measured [33] the 
single diffractive cross section: ctsd = 9.46 ± 0.44 mbarn. The three quantities cr to t, °<=i ano - ^ci are related by the 
unitarity relation (|39[) . and therefore only two of them are independent. In the following we will fix our attention on 
atot and B ci . 

Our formalism allows only the calculation of the total diffractive cross section, summing over single and double 
diffraction processes. In order to compare the model to the single diffraction measurement of CDF wc have therefore 
to include some estimate of double diffraction. We will use the result ((76)) that allows to estimate the complete 
diffractive cross section Odiff from the measurements of the elastic and single diffractive one as: 



This hypothesis leads us to estimate a^is at y/s = 1800 GeV from the CDF data as approximately 11.6 mbarn. 

Neglecting the measurement of the diffractive cross section and considering only the experimental results for cr to t 
and -B e i, there is an infinity of sets of parameters (cr c ik, r , w ) that reproduce the data. This infinity of solutions can 
be parametrized by the value of w, that can take any non-negative value. 

The limiting case w = corresponds to the simple eikonal model. For w = the values of the parameters that 
reproduce the central value of the CDF measurements are: er C ik — 124.1 mbarn and ro ~ 0.2527 fm (~ 1.08 R p ). To 
this (w = 0) solution corresponds a vanishing diffractive cross section. 

The solutions [w, <7 e ik(w), ro(w)] that reproduce the central value of the CDF measurements for crtot and B e \ at 
y/s = 1800 GcV are shown in fig. [3J Increasing w the value of <7 e ik(w) grows monotonically, while ro(u>) decreases. 
The triplets {w,a e i^(w),ro(w)} result in identical crtot and B c \, but produce different diffractive cross sections, with 
fdiff(w) growing monotonically with w as also shown in fig. [21 The value Cdiff = 11.6 mbarn is obtained for w = 3.48. 

Perhaps the most striking feature of figure[3]is the rapid increase of cr c ;k with w. As an illustration, for w ~ 3 (that 
results in Odiff — 10.7 mbarn) one needs a smaller ro (ro = 0.186 fm) and cr e ik — 580 mbarn, that is almost 5 times 
larger than the value of a e ik that reproduces the measurements for w = 0. 

For a qualitative understanding of these results it can be instructive to consider figure SJ Curve (a) shows the profile 
function r e i(6) that corresponds to the solution with w = that we have just discussed. Curve (a') shows the profile 
that is obtained for the same parameters (tr e ik, fo), that is for the same (n(6, s)), of the w — solution, but using the 
value w = 3. The resulting profile function is smaller (that is produces a smaller crtot) and broader (implying a larger 
-B e i). These features can be readily understood from inspection of equation ([ST]) . In order to obtain the desired values 
of <7tot and B c \ using the model (|8ip and w = 3 one needs to modify the quantity (n(b, s)), choosing both a larger 
cr c ik to increase the area under the profile, and a smaller ro to obtain the desired value of (b 2 ) oc B. The solution is 
shown in fig. 2] as curve (b). The profile functions of curves (a) and (6) in fig. [1] produce identical crtot, and identical 
da e i/dt for small \t\, but differ in the description of elastic scattering at large \t\. 

To summarize this discussion: in the simple eikonal model the impact parameter multiplicity distributions (n(b, s)) 
uniquely defines the profile function. In our model the profile function is determined (see equation (|8ip) also by the 
function p(a), and different choices for the shape of p(a) produce different profiles, and therefore different values of 
the total and elastic cross sections. Viceversa, the estimate of (n(b, s)) (or in terms of the parametrization ([43)) the 
values of cr ik and ro) that reproduces the measured values of crtot and cr G i (or Otot and -B i) strongly depends on the 
assumptions made for p(a). 

The function p(a), and in particular its width, controls the size of the inelastic diffractive cross section, therefore 
one can obtain information about its properties from the experimental data on the rate of diffractive events. The 
bottom line is that it is essential to include in a consistent way inelastic diffraction in the theoretical framework that 
describes hadronic cross sections. 

These considerations are the main qualitative results of this work: the consistent introduction of inelastic diffraction 
in the eikonal formalism results in: 

1 . an eikonal cross section a e ik that is several times larger than estimates based on the simple' eikonal model that 
does not consider explicitly diffraction; 

2. a narrower distribution of hadronic matter. In the case of protons, this distribution is estimated as narrower 
than the charge distribution infered by the electromagnetic form factor. 




(85) 
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These effects can have important consequences in the prediction of the properties of particle production in high 
energy hadron interactions, if one takes into account the interpretation of the eikonal as a description of the multiple 
interaction structure of the collision. In this case the ratio er e ik/ 'omei has the physical meaning of the average number 
of elementary interactions per inelastic event, therefore our results imply that this average number of elementary 
interactions is several times larger than previous estimates. The precise way to relate this quantity (the average 
number of elementary interactions per collision) to observable quantities, such as the multiplicity distribution, depends 
on a number of additional assumptions that have to be made in a Montccarlo modeling of multiparticlc production. 

The theoretical framework we are considering predicts not only the average number of interactions in a collision, 
but also the detailed multiplicity distribution for such interactions in one collision. 

In section [TT] we have given in equation (|26p the multiplicity distribution of the number of hard observable jets per 
collision, in terms of the quantities (nj et (6, s)) and p(a). The generalization to the multiplicity distribution of the 
total number of elementary interaction in a collision can be immediately obtained replacing the average number of 
hard interactions at a fixed impact parameter and cm. energy (rij e t{b,s)) with the average for the total number of 
elementary interactions (n(b,s)). In addition, the most economical assumption is to assume that the functions p{a) 
relevant in the two cases are (at least approximately) equal. 

For the simple functional form of p(u) given in equation (|77[) the integrals over a in the analogous of equation Q26|) 
can be performed, with the result: 

J d 2 b (n(b, s)) k {1 + w (n(b, s))r( fe+ ») (86) 

This distribution of the number of elementary interactions should be inserted in Montecarlo implementations to have 
predictions for the charged particles multiplicity and other observable quantities. 

VI. ENERGY DEPENDENCE 

The model we have outlined in the previous sections consider the relation between directly observable quantities 
such as the total and elastic cross sections and on the other hand the eikonal cross section cr e ik and the distribution 
of hadronic matter in the colliding particles (simply parametrized by the quantity ro), with an additional parameter 
w that is related to the width of the fluctuations in the parton configurations of the colliding hadrons. Taking into 
account these fluctuations allows a consistent treatment of inelastic diffraction. The values of <7 e ik and ro that correctly 
describe the data, have a strong dependence on the parameter w, and therefore on the measured values of the cross 
section for inelastic diffraction. 

A calculation of the evolution with energy of the hadronic cross section requires additional theoretical assumptions 
to predict the energy dependence of the model parameters. To gain insight on this problem, we have taken a 
phcnomcnological approach and we have considered a representative subset of the available high energy data. A few 
high energy experiments have measured both the total and elastic cross section, together with the forward slope B e i- 
These results can be described in terms of our 3-paramcter {<7 C ik, ro, w} model, using the same approach discussed 
for the CDF data at yfs = 1800 GeV. The results are shown in fig. [5] and [51 In these figures the points give the values 
of ro and <7 c ik that reproduce the measurements of the pairs (otot, B) using two assumptions for the third parameter: 
w = and w = 3. The errors on the estimates reflect only the experimental statistical errors. At yfs = 1800 GeV, 
one has two independent measurements of the total cross section by the CDF [2(| and E710 [38[ experiments, both at 
the Fermilab pp collider. The point at yfs = 546 GeV is also from CDF, while the point at s/s — 62.3 was obtained 
for pp scattering at the CERN ISR collider 

The calculation with w = corresponds to the simple eikonal model, and in the framework to our model leads to a 
vanishing inelastic diffraction cross sections. For each pair of experimental results ((7tot,-B) we have performed a scan 
of w similar to the one that we have shown in detail for the CDF point at s/s = 1800 GeV (see also fig. [3]), calculating 
the pair (ro(w),a e a i (w)) that reproduces the experimental results. This also imply a value <Jdm(w) obtained from 
equation ([84]) . 

At yfs = 1800 GeV both the CDF [37} and E710 [38[ have measured the single diffractive cross sections. Averaging 
with equal weight (to take into account large systematic uncertainties) and using the ansatz (|76j) to estimate double 
diffraction we have estimated a dm — 10.7 mbarn. This value of the diffractive cross section is reproduced in our model 
(at the corresponding energy y/s = 1800 GeV) with w ~ 3. Including a 20% uncertainity on the estimate of Cdiff , w 
can be estimated at this energy asw = 3+J;g. 

The comparison of the calculated diffractive cross section with the data is problematic because the discrepancies 
between the different experimental results (see fig. [7]) clearly indicate the presence of significant systematic errors, 
moreover, as we have discussed before, one has the theoretical uncertainty related to the ratio between the single and 
double diffraction contributions. 
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Our numerical studies indicate that the choice of an energy independent value w ~ 3, that reproduces the diffractive 
cross section measured at the Fermilab collider at = 1800 GeV, gives in fact a reasonably good description of the 
experimental results on diffraction at all energies. The assumption of an energy independent value of w is clearly the 
simplest one, and in view of the fact that it produces a reasonable agreement with the available data it will be made 
in the following. 

It is interesting to note (see fig. [5]), that assuming for w the constant value w = 3 the resulting values of the 
parameter r$ are also consistent with an energy independent value r$ ~ 0.19 fm. In the simple eikonal model, as 
discussed before, in order to reconcile the growth of Ctot and B c \, it is necessary to increase the width of the overlap 
function A(b, s) with s, and therefore (using our paramctrization) to increase ro(s). The necessity of this growth is 
evident in fig. [5] 

The model we are discussing requires an overlap function that is first of all significantly narrower than previous 
estimates based on the simple eikonal model; moreover (and in contrast to the simple eikonal model) the overlap 
function is energy independent. It may appear surprising that the overlap function is narrower than what is estimated 
on the basis of the proton charge distribution. A possible explanation is that the dominant contribution to the overlap 
function is the scattering between soft gluons. The narrow A(b) predicted in this model therefore implies that the 
impact parameter distribution of soft gluons is (i) narrower that the charge distribution (that is presumably controlcd 
by valence quarks), (ii) independent (for small x) from the x of the gluons. It should soon be possible to test these 
hypothesis with studies of the impact parameter PDF's. 

Figure [S] shows the energy dependence of the third parameter of our model, <J ik(s). There are two remarkable 
features in this behaviour. The first is that (as we have already discussed in section IVA"]) the values of cr e ik needed to 
describe the experimental data in a model that includes diffraction are significantly larger than estimates based on 
the simple eikonal model, the second is that the growth of cr e ik(s) with cm. energy is significantly more rapid. 

The two straight lines in fig. [5] are power law fits of the estimated values of form K s a . For the simple eikonal model 
(w = 0) the power law fit is er c ik(s) = 23.8 s 0A0 mbarn (with s measured in GeV 2 ), for the best fit model (w = 3) the 
power law fit is er C ik(s) = 29.7 s 018 mbarn. 

The choice of a power law fit is however clearly not necessary, and the extrapolation of the fit at higher (and lower) 
energy is therefore very uncertain. Motivated by fits of the er e ik(s) = "soft + 0hard(s) that include an (approximately) 
constant soft component and an energy varying hard component, we have fitted the data with the form <r ik(s) = 
o"o + K s 35 , obtaining the result <7 e ik(s) = 95 + 2.1 s 0,35 mbarn (and s measured in GeV 2 ). 

The motivation for the functional form of this fit (that should however also be considered as purely phcnomcnolog- 
ical) is that integrating above an energy independent p" lln , the jet cross section crj et (p™ m , y/s) (see equation ([3])) has 
qualitatively the behaviour: 

^jctbr\ V^) OC j-^ L_JW a s e log s (87) 

where t = 4 (p" lm ) 2 /s, and the quantity e is related to the behaviour of the PDF's for x — > 0: 

lim /(x) ~ (88) 

x— >0 X ^ 

The behaviour i[8?]) can be easily be obtained from the convolution of PDF's with the asymptotic form (|88|) . Recent 
measurements of the PDF's at HERA [39|, H3, 51], |42j have shown that their behaviour for x — ► can be reasonably 
well represented with the functional form (|8"8")) and e ~ 0.3. 

It is interesting to note that it has been argued that the very fast growth of the jet-cross section with y/s implied 
by the small x behaviour of the PDF's is problematic, and in fact unphysical. Note that in the simple eikonal model, 
an energy dependence of c e ik(s) of type s 3 (or faster) is not acceptable because it implies that CTtot(s) grows with 
energy more rapidly than the observations. Following this observation, the rapid growth of <7j et (s) with s has been 
tamed assuming that the threshold p" lm of applicability of perturbation theory also grows with s. This growth has 
been connected to phenomena of "saturation" , or screening among the partons. In the framework of the model we 
are considering a fast growth of cr ik(s) is not only acceptable but in fact necessary. A simple model where the growth 
of cr e ik(s) is explained with the dominant contribution of a minijet cross section, calculated perturbatively above an 
energy independent p" lln can provide a tr e ik(s) with the needed properties. For consistency, it is however necessary 
that the effects of saturation and parton screening are small. 

The model we are describing, already at y/s = 1800 GeV has a ratio Ueik/cinei — 10. This implies that the number 
of elementary interactions in an inelastic collision at this energy is also approximately 10. This, at first sight, may 
appear too large. The potential danger is that this large average number of elementary interaction per collision could 
result in a too large average multiplicity and in a too soft inclusive spectrum of particles in the final state. These 
questions can (and should) be addressed properly with a detailed Montecarlo calculation, that includes a modeling of 
particle production in the presence of different numbers of elementary interactions. 
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Figure [7] shows our calculation of the diffractive cross section Odiff(s) including the extrapolation to high energy. 
The calculation is performed with equation (|84[) . using energy independent values w = 3, r = 0.19 fm and the two 
paramctrization of cr G ik(s) shown in fig. El 

Figure[5]shows the result of our model for the total cross section, comparing with the available data and extrapolating 
at higher energy. In the figures we show two calculations based on equation (|82[) using (as in the previous figure) the 
constant values w = 3 and tq = 0.19 fm, and the two parametrizations of cr e ik(s) (the results for the model are only 
plotted for > 60 GcV). 

Figure [8] also shows the parametrizations for cr^ t (s) and af^ t (s) suggested in the PDG [ID, [44| . The PDG estimate 
of the extrapolation of the total pp cross sections falls in between our two estimates, that mark a range of uncertainty 
in our prediction. 

At the LHC energy (y/s = 14 TcV) the PDG prediction is ertot = 112.2 mbarn while our two calculations give 
Ctot = 98.1 and 120.8 mbarn [at ^fs = 10 TeV the PDG predictions is otot = 105.7 mbarn, our calculations give 94.2 
and 112.7 mbarn]. The extrapolation to y/s = 4.33 x 10 5 GeV (that corresponds to a proton cosmic ray particle with 
energy = 10 20 eV is ertot = 194 mbarn for the PDG extrapolations, and 146 and 229 mbarn for our calculations. 

Recently two groups [H, Q have discussed predictions of the total pp cross section at high energy in the framework 
of models that include a treatment of diffraction with a multi-channel cikonal and multi-pomeron interactions. Both 
groups arrive to a similar conclusion, namely that the inclusion of diffraction reduces the estimate of the total cross 
section at LHC energy. For the two groups the estimate of er to t at yfs — 14 TeV is of order 90 mbarn, approximately 
20% smaller than the PDG. Their estimates of the total cross section grows very slowly with energy reaching ertot ~ 
108 mbarn (for '9\) and cr tot ^ 98 mbarn (for 0) at s/s = 10 5 GcV (that corresponds to E uh = 5.33 x 10 18 eV). 

The main point that we want to make here, is that it is certainly the case that given a model for (n(b, s)} (that 
in the simple eikonal model is simply equal to twice the eikonal function \(b, s)), the inclusion of diffraction reduces 
the cross section. However, the estimate of the total cross section and of its dependence on energy also involves the 
calculation of the function (n(b, s)). 

The conclusion that the cross section at LHC is of order 90 mbarn obtained by the authors in || Q should not be 
considered as a consequence of the inclusion of diffraction in the theoretical framework, but rather as the consequence 
of the entire set of theoretical assumptions of their models. 

Figure [5] shows, plotted as a function of -y/i, the predictions of our model for the slopes B e \(s) and i?diff( s ) of the 
differential cross sections do~ e \/dt and dadis/dt. For each slope, the figure shows two curves that differ for the use of 
the two different parametrizations of (T e ik(s) that are shown in fig.0 and have already been used in fig. [7] and [5J Note 
how -Bdiff(s) is always larger than B e \(s). 

VII. SUMMARY AND OUTLOOK 

In this work we have discussed the problem of multiple parton interactions in hadron collisions. If one takes into 
consideration only the parton scatterings that have sufficiently large momentum transfer, it becomes possible to detect 
the final state partons as high p± jets, and determine event by event the number of hard parton interactions that 
are present. It becomes therefore possible to study experimentally the multiplicity distribution of parton interactions 
above for example p™ ln . At sufficiently high cm. energy the probability of having more than one high pj_ parton 
scattering in a single collision can be appreciable, and the hard scattering multiplicity distribution becomes non trivial. 

The calculation for the inclusive distribution of high p± parton scattering is a textbook application of perturbative 
QCD and can be performed from a knowledge of the standard PDF's. A theoretical prediction of the multiplicity 
distribution of high p± interactions in high energy hadron collisions is however a highly non trivial problem that 
requires the introduction of additional theoretical concepts. 

The standard PDF's give the inclusive probability density for finding one single parton with fractional longitudinal 
momentum x. This probability density is obtained integrating over the parton transverse momentum, and integrating 
over all possible momenta of the other partons in the hadron. If one wants to compute the probability to have exactly 
n hard interactions in one collision the information contained in the PDF's is clearly insufficient. One needs to know: 

1. the probability to find a parton of a given x at different impact parameters with respect to the hadron center 
of mass; 

2. the correlated probabilities for finding different partons at {x\, b\), (x2, ^2), (^3, ^3) , • ■ ■• 

The first problem should be addressed introducing impact parameter dependent PDF's, F^x, b, Q 2 ) that give the 

probability of finding the parton of type j with fractional longitudinal momentum x and impact parameter b probing 
hadron h at the scale Q 2 . This problem has not yet a well determined solution, and all studies of this problem have 
made the simplification to assume that the dependences on x and b of the impact parameter PDF's factorize, that 
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is: Ff(x,b,Q 2 ) = fAx 




Q 2 ) p(b)i and estimated p(b) with simple phcnomenological considerations. Studies of the 
131 fl4| should soon be able to shed light on this question. 



Generalized PDF's h 



The problem of obtaining correlated PDF's that give the probability to find simultaneously several partons in 
different elements of phase space is clearly much more difficult and complex. In this work we have suggested to 
parametrize the effects of our lack of knowledge about the correlated PDF's introducing the "effective configuration 
probability distribution" , that is one function p(a) the real, positive variable a. Each one of the configurations of 
partons in the pair of colliding hadron has associated the real number a. The physical meaning of a is that the 
expected number of parton interactions that corresponds to the parton configuration C is n(C) = (n) a, where (n) is 
the average over all configurations. The first two moments of the function p{a) are unity (because of the normalization 
of a probability density and to reproduce the correct (n)); increasing the 2nd moment of the p{a) distribution the 
width of the multiplicity of parton interactions grows. 

If one considers not only a subset of detectable (high p±) parton interactions, but all of them, it becomes possible to 
relate these elementary interactions with the total and elastic cross sections. This general idea has been implemented 
in many works using an eikonal formalism. A crucial ingredient of these models is the quantity (n(b, s)) that gives 
the average number of elementary interactions for a hadron collision at impact parameter b and cm. energy ^/s. In 
the simple eikonal model this quantity is related to the clastic scattering profile function by the relation r e i(6, s) = 
1 — exp[— (n(b, s))/2]. The corresponding inelastic cross section is then: 



The physical interpretation is that an inelastic interaction corresponds to absorption and to at least one elementary 
interaction, assuming Poisson fluctuations in their multiplcity. The same considerations that wc have outlined for 
hard interactions however apply, and it is natural to expect that fluctuations in the number of elementary interactions 
n are in fact much broader than poissonian because of fluctuations in the "configurations" of the colliding hadrons. 
This effect can again be parametrized with a function p(a). For example the inelastic cross section can be rewritten 
as: 



One can see that the parameter a controls the "transparency" of a hadron collision. Different "configurations" of 
the colliding hadrons have transparencies that are related to a. Good and Walker [U have proposed that inelastic 
diffraction originates from the different absorption of the different components of the colliding hadrons. Therefore our 
formalism can be applied to the calculation of the inelastic diffractive cross section, and in fact unavoidably implies 
the presence of inelastic diffractive processes. 

In other words, the function p(a) allows to relate the quantity (n(b, s)) to the total and elastic cross sections, and 
at the same time fixes the value of the diffractive cross section. Viceversa, from the data on the total and elastic 
cross section, together with the data on inelastic diffraction it is possible to extract information on (n(b, s)) and on 
the properties of p(a). 

We have performed an analysis of the data on pp and pp collisions obtained at high energy colliders, and obtained 
information on (n(&, s)) and p(a). For the study of the properties of p(a) we have used as a first approximation a 
simple analytic form that depends on a single parameter. 

To describe the measured diffractive cross section one is forced to have a fuction p(a) with a large variance. This 
in turn has very important consequence on the parameters that describe (n(b, s)). It is remarkable that wc find that 
(within significant uncertainties) the function p(a) is independent from energy; moreover parametrizing (n(b, s)} in 
the form: (n(b, s)) = cr ik(s) Mp, s) as the product of an eikonal cross section times a geometrical overlap function, wc 
find that the geometrical factor can be taken as energy independent, in contrast with results obtained in the simple 
eikonal model that neglects fluctuations. The eikonal cross section <7 e ik(s) is much larger and grows much faster with 
energy than in the simple eikonal model. Such a rapid growth can however be readily explainable assuming that it is 
controled by the increase of (Tj e t(p™ m , s) with s assuming a constant p™ m and negligible screening effects. 

We note that, at least in first approximation, the function p{a) that we have extracted from the study of the total, 
clastic and diffraction cross sections, is also applicable to the study of the multiplicity distribution of high p± jets. 
One can therefore make the prediction that the distribution of the number of hard interaction per event will be broad, 
with a non negligible number of events containing several interactions. 

The prediction of the total cross section at LHC depends on the energy dependence of tr e ik(s). This is a problem 
we have not discussed in detail here. It seems however natural to expect a result around 110 mbarn with however 
a significant uncertainty. In the model we are discussing however the eikonal cross section cr e ik(s) is large and since 
Ceik(s)/<7i n oi(s) is equal to the number of elementary interactions in a collision one is lead to expect a large charged 





da p(a) {1 — exp[— (n(b, s)) a}} 
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particle multiplicity and a soft inclusive spectrum in the final state. These consideration are also relevant for the 
study of ultra high energy cosmic ray showers. 
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APPENDIX A: ELECTROMAGNETIC FORM FACTORS AND OVERLAP FUNCTION A(b) 

In [l[ Durand and Pi estimated the overlap function A(b) for pp collisions from the electromagnetic form factor of 
the proton. The simple physical idea behind their derivation is that the overlap function is the energy independent 
geometric overlap of the hadronic matter distributions in the colliding particles. More cxplicitcly, one defines the 
spatial distribution of matter in the hadron h as Ph{f), with the normalization condition: 

d 3 r p h (r) = 1 (Al) 

The density in the transverse plane is then be obtained with a simple integration: 

Ph(b) = f + dz Ph (vV + z 2 ) (A2) 

J — OO 

The overlap function in the collision between hadrons hi and /12 is then obtained as: 

A(b) = f d 2 b x d 2 b 2 p h ,{bx) p h2 {b 2 ) 8[b-h + h] (A3) 



J d 2 b A{b) = 1 (A4) 



The normalization condition 



follows automatically from the normalization of ph x {r) and ph x (r) given by (jAlj) 

To estimate the density of Phif) of hadron h, Durand and Pi make the assumption that it is simply the Fourier 
transform of its electromagnetic form factor. For the proton one has: 

W) = ( TT i f?F (A5) 

with R p = 0.234 fm (or R~ 2 = 0.71 GcV 2 ) and correspondingly p p (r) cx e~ r ' Rt '. The geometric convolution for 
proton-proton collisions is: 

A(b) = &3 z K 3 ( —\ (A6) 



where Kz(x) is the modified Bcsscl function of the third kind. 



APPENDIX B: FOUR CHANNELS MODEL 

It is instructive to discuss a "minimum" model that implements the Good- Walker ansatz 0] for inelastic diffraction 
in the collision between two hadrons, where all calculations can be performed cxplicitcly The minimum model has 
4 channels, to describe the 4 possible types of scattering (target, projectile and double diffraction together with 
elastic scattering). Without loss of generality we can consider the scattering irp (here 'V and "p" are labels to 
represent arbitrary hadrons). Each of the two colliding hadrons can undergo inelastic diffraction with a transition 
to an additional state. We will label the excited states for the projectile and target particles as u p" and "A"; one 
therefore has to consider the transitions: n — > tt* = p and p — > p* = A. In the 4-channel model one has to study the 
4-dimcnsional vector space spanned by the orthonormal basis of the 4 physical states \ip m )'- 



{|^m)}(m=i,4) = {\np),\nA),\pp),pA)} 



(Bl) 
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One can (in principle) study the 4x4 transitions ((pf\S\tfi}. In practice of course one is limited to the study of 
the transitions \irp) — ► \tp m ) (that correspond to the processes of elastic scattering, target, projectile and double 
diffraction). The scattering amplitude is a 4 x 4 matrix: 



Z 7T 



The differential cross section for the transition i — > / is 



dt 



(t,s) = TT 



HZ*) 



fi 



Integrating over all t values one obtains the transition cross sections 



(B2) 



(B3) 



T(b,s) 



n 



The profile matrix T can be written in terms of the eikonal matrix x(6, s): 

t(b, s) = l- cxp [-x{b, s)} . 
Using the Good and Walker ansatz, the eikonal matrix x(b, s) takes the form: 



X(b, s) 



i(b,s)) 



M 



(B4) 



(B5) 



(B6) 



where (n(b, s)) has the usual meaning of the average number of parton interactions in a np collision, and we have 
introduced the 4x4 matrix M. This matrix must be real and have 4 real and positive eigenvalues ctj, moreover one 
must have M\\ = 1. 

This is a consequence of the fact that one can define (as in the previous section) the states \ipj) as the eigenstates 
of the M matrix. These states undergo only absorption or elastic scattering, and each has a "transparency" Pq = 
e ~{n(b,s)) oij ^ therefore (n(6, s)) ctj can be interpreted as the average number of interactions for the state \ipj) and 
therefore this quantity (and Oj) must be positive. Moreover, to have the correct average multiplicity of elementary 
interactions for the intiali state \tpi) = \ttp) one must have: 



= i 



(B7) 



that implies M\^\ = 1. 

The matrix M can be constructed cxplicitcly making the additional hypothesis that the (4-dimensional) space of 
the physical states is the direct product of two (2-dimensional) spaces for the beam and target particle, and moreover 
that one has time reversal symmetry, and the amplitude for the transitions ir — > p (p — > A) and p — > tt (A — * p) are 
equal. With this assumptions the most general form for the matrix M is: 



M 



1 Ar 



1 P P 
d v l-2e ? 



1 - 2e ff 0„P P p p (l-2e 7r ) 

P P P.Pp 1 - 2e p p 7v (l-2e p ) 

\P„P P P P (1 - 2e n ) & (1 - 2e p ) (1 - 2e OT ) (1 - 2e p ), 



(B8) 



The eigenvalues and eigenvectors of the matrix M are easily calculable, noting that each of the 2x2 matrices of 
form 

1 P 
P l-2e 
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has eigenvalues: 



Ai,2 = 1 ± 7 



where 



and the corresponding eigenvectors are: 



7= 



vi, 2 = -j= {±vT±r, y/l Tr} 



(B9) 



(BIO) 



(Bll) 



e/7 = 



with ' 

The eigenvalues of the 4x4 matrix M are then: 

ay = (1 ± - e w ) (1 ± j p - e p ) 
(with j G {1,2,3,4}). The condition that the eigenvalues are non-negative gives: 



Eff , p < 1/2; Pl p 



< l-2e 



7T,p 



(B12) 



(B13) 



The rotation matrix C rn j that connects the scattering cigenstatcs \ipj) to the physical states \ip m ) is: 



'>i j 



( yTTT^ yi~r7~; a/T+~ 



y i + a/ i - ?> 



V I + ry a / 1 + ?> 



V I + r n y / l - r p y /l -r^ y /l- r p \ 
— y /l - ?V a /1 - ?> a/ 1 + ry ^ /l - r , 

\ZTT7~ v /TTr7 y/T~T^y/TT- 



(B14) 



rp -VI + r n y/l - r p -y/l - ry y/l + r p y/l +r IA /l + r p ) 
The profile functions for the different scattering processes can now be given explicitly as: 

r -Vc -c fi ex P f 



(B15) 



The model outlined above requires an estimate of the function (n(b, s)) that can be interpreted as the average number 
of "elementary" interactions for a hadron crossing at impact parameter b and cm. energy y/s. In the general case of 
the collision of two different hadrons (such as in ir p scattering) the model has 4 additional parameters e ni (3 P , e p ) 
that describe the matrix structure of the cikonal function. Obviously for pp scattering the model has only two 
parameters (/3 P , e p ). 

With the labeling of the physical states that we have been using (namely: \ipi) = \np), \if2) = |vrA), l^) = \pp) 
and 1 1/34) = |pA)) the integration over all impact parameters 6 of |r 1:L (rb)| 2 yields the elastic cross section, the integral 
of |r2i(6)| 2 (|r3i(6)| 2 ) gives the target (projectile) single diffraction cross section, and finally the integral of |r4i(6)| 2 
gives the double diffractive cross section. 

To connect this analysis to the discussion performed in the main text, we note that we can define the 4 quantities 
Pj that are the probabilities Kv'ilV'j')! 2 to find the initial state |y>i) = \ttp) in the scattering eigenstates \ipj). The pj 
are given by: 

Pj = KV'jKp)! 2 
It is straightforward to verify that: 



( lv ±e ) (l P ±e P ) = l 
4 7 T 7 P 4 
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One can now define the function p(ct): 



= 5[a-a.j\ 



(B18) 
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This function, as a consequence of equations (|B17[) satisfies the conditions: 

/>oo />oo 

/ da p(a) = 1, / da a p(a) = 1. 

Jo Jo 

It is now straightforward to see that one can recast the expressions for the total, elastic, absorption and diffractive 
(that is the sum of the target, projectile and double diffraction) cross section as integrals over a identical to the 
expressions (|55H52^) in section HV Bl 

It can be interesting to note that the 2nd moment of the p{a) distribution is given by: 

J da a 2 p(a) = J^pj a 2 = (1 + p p ) (1 + (B19) 
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FIG. 1: The points are measurements of the total cross section cr to t and of the forward slope B B \ of the elastic scattering 
for pp and pp collisions at collider energies. The lines correspond to predictions based on the simple eikonal model using the 
parametrization of equation (|45[) for (n(b,s)). The three lines are computed for three values of the ro parameter (ro = R p , 
1.1 R p and 1.2 R p ). The ISR data at y/s = 52.8 and 62.3 Gev is from %2$; the UA1 data at y/s = 540 GeV from Hi|; the CDF 
data at = 546 and 1800 GeV from [H,^; the E811 data at 1800 GeV from 
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FIG. 2: Plot of the function p(a) given in equation (J77J) for four values of the parameter w (w — 0.05, 0.5, 1, 3). For w — > 
the function takes the form 8 [a — 1]. 
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FIG. 3: The middle and bottom panel show the triplet of parameters (w, a c ^(w) and ro(w)) that re pro duce (using equations 
dH2J an d (H31) with expression (gSJ for (n(b, s))) the measurements of <7 to t and B cl obtained by CDF [H, [23] at ^/s = 1.8 TeV 
(note the logarithmic scale in the bottom panel for cr ik). The top panel shows the corresponding value of the diffractive cross 
section. 
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FIG. 4: Profile function r c i(6) for pp scattering as a function of the impact parameter b. The curve (a) is calculated in the 
simple eikonal of equation (|4ip . using for the {n(b, s)) the parametrization (|45[) with a c ;k = 124 mbarn and ro = 0.253 fm. The 
corresponding values of <7 t ot and B e \ are a to t = 80.3 mbarn and B c \ = 16.98 GeV~ 2 . Curve (a) is calculated assuming the 
same interaction profile (n(b,s)) (that is the same parameters and ro) as for curve (a) but using the model of equation 
(|81|) for the profile with the form (|77[) for p(a) with w = 3. The resulting profile function is smaller (implying a smaller Ctot) 
and broader (implying a larger B e \). The profile (6) is calculated using the same model used for curve (a') with the same value 
w = 3, however the parameters that describe {n(b, s)) are now <r c ik = 582 mbarn and ro = 0.186 fm. The profile function (6) 
results in the same atot and B c \ as curve (a). 
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FIG. 5: Values of the ro parameter that reproduce the experimental data for otot and B c i obtained at the ISR pp collider 
(Vs = 62.3 GeV), and at the Tevatron pp collider (y/s = 546 GeV by CDF, and i/a = 1800 GeV by CDF and E710). The solid 
(empty) points are calculated for w = (w = 3). The dotted line at ro = 0.234 fm corresponds to the proton charge radius R p . 
The dashed line corresponds to the constant value ro = 0.19 fm, and is a reasonable representation of the results for w = 3. 
The corresponding values of cr c ik are shown in fig. [6] 
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FIG. 6: Values of the cr e ik parameter that reproduce the experimental data for a to t and B c \ obtained at the ISR pp collider 
(Vs = 62.3 GeV), and at the Tevatron pp collider (y/s = 546 GeV by CDF, and y/s = 1800 GeV by CDF and E710). The 
solid (empty) points are calculated for w = (w = 3). The corresponding values of ro are shown in fig. [5] The dashed line is 
a power law fit (a e ik(s) = K s a ) to the results for w — 0. The thin (black) line is a fit to the results for w = 3 with the same 
power law form. The thick (blue) line is a fit to the same points with the form a e ik(s) = <Jo + K s a with a — 0.35. 
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FIG. 7: Inelastic diffraction cross section calculated according to equation (|84l) using constant values w = 3 and ro = 0.19 fm. 
For the thin (black) [thick (blue)] curves we have used for cr e ik(s) the fit shown with the corresponding lines in fig. [HI The 
experimental results are for single diffraction only (Schamberger jig, Armitage [Hj], UA4 [I?]], UA5 CDF [371, E71 ° H3l)- 
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FIG. 8: The points are measurements of the pp and pp total cross sections. The dashed lines are the fit of a to t{s) suggested in 
the PDG [3]. The other two lines are predictions obtained from equation (|82|l using constant values w = 3 and ro — 0.19 fm. 
For the thin (black) [thick (blue)] curve we have used for <r c ik(s) the fit shown with the corresponding lines in fig. [6] 
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FIG. 9: Slope at t = of the differential cross sections da e \/dt and dadis/dt for elastic and inelastic diffractive events. 
The predictions are calculated with equations (|65|l and (|66[) . using the functional form (|77|l for p(a) with ui = 3, and the 
parametrization (|45f) for (n(6, s)} with ro = 0.19 fm. The solid and dashed curves use the two parametrizations of cr e ik(s) shown 
in fig. El 
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